Preprint typeset in JHBP style - PAPER VERSION 



ZU-TH 04/13, LPN13-012 



Two-loop QCD helicity amplitudes for gg ^ Z g and 

9 9 Zj 



Thomas Gehrmann, Lorenzo Tancredi, Erich Weihs 

Institut fiir Theoretische Physik, Universitdt Zurich, Wintherturerstrasse 190, 
CH-8057 Zurich, Switzerland 



Abstract: We compute the helicity amplitudes for the processes gg — )• Zg and gg — )• Z7 
to two loops in massless QCD. The perturbative expansion of these processes starts only 
at the one-loop level, such that our results are a crucial ingredient to the NLO corrections 
to Z'j and Z+jet production through gluon fusion. 

Keywords: QCD, Colhder Physics, NLO and NNLO Calculations. 



Contents 



1. Introduction 1 

2. Kinematics and notations 2 

3. The general tensor structure 5 

3.1 The gauge-fixed tensor structure 7 

4. Helicity amplitudes 9 

4.1 gggV: The amplitude in spinor helicity notation 9 

4.2 ggjV: The amplitude in spinor helicity notation 10 

4.3 Analytic continuation to the scattering kinematics 11 

5. Outline of the calculation 13 

5.1 UV Renormalisation and IR subtraction 14 

5.2 Simplification using the Symbol formalism 15 

6. Checks on the result 16 

7. Conclusions and Outlook 17 

A. One-loop helicity amplitudes 18 

A.l V ^ ggg at one loop 18 

A. 2 V — )• gg^ at one loop 18 

B. Two-loop amplitudes: all-plus helicity coefficients 19 

B. l V ggg at two loops 19 
B.2 V — )• gg^ at two loops 23 



1. Introduction 

The production of vector bosons at hadron colliders is, to a first approximation, induced by 
quark-antiquark annihilation. Including corrections from higher orders in the perturbative 
expansion in QCD, other processes will also contribute to vector boson final states. These 
contributions are suppressed by higher orders in the strong coupling constant ag, but 
could receive a numerical enhancement through the relevant parton-parton luminosity. In 
particular, in high-energy proton-proton collisions at the LHC, gluon- induced higher-order 
processes can become of comparable importance to quark-induced processes due to the 
large gluon luminosity at invariant masses relevant to vector boson production. 
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Vector-boson production in gluon-gluon collisions is mediated through a quark loop, 
which vanishes for the exclusive gg ^ V vertex due to Furry's theorem. The gluon-gluon- 
induced subprocess becomes relevant for the production of vector boson pairs (VFVF, ZZ, 
77 and Z7), or for the production of a neutral vector boson and a gluon. The leading- 
order scattering amplitudes for these processes all involve a closed quark loop. The result- 
ing gluon-induced contributions from one-loop squared [1-3] processes (that appear only 
at next-to-ncxt-to-leading order in the formal perturbative expansion of the full process) 
were evaluated a long time ago [4-8], and typically found to yield a contribution that 
amounts to 10-20% of the total cross section. Inclusion of these gluon-gluon subprocess 
contributions often results in an enhanced theoretical uncertainty on the prediction, since 
the one-loop squared process is effectively Born-level for this combination of partons. To 
stabilise these predictions, the computation of the next perturbative order in vector-boson 
pair production or vector-boson-plus-jet production in gluon fusion is required. Techni- 
cally, such a calculation amounts to computing the corrections from single real radiation or 
single virtual exchange to the Born processes. With the Born process itself being a one-loop 
amplitude, one thus requires the two-loop corrections to the relevant partonic amplitudes. 
Up to now, these were obtained [9] only for gg — )■ 77, where the NLO correction to the 
gluon-induced process was found to be sizeable and important in the stabilisation of the 
theoretical prediction for photon pair production [10,11]. 

In this paper, we derive in massless QCD the two-loop corrections to the helicity 
amplitudes relevant to the production of a Z-boson in association with either a real photon 
or a hadronic jet in gluon-gluon collisions: gg Zj and gg — > Zg. For these processes, the 
one-loop amplitudes involving an extra gluon in the final state can be obtained using by- 
now standard methods for the computation of one-loop multi-leg processes [12-17]. With 
the results derived here, a complete NLO calculation of Z7 and Zj production in gluon 
fusion becomes thus feasible. 

This paper is structured as follows: in Section 2, we fix the notation and discuss the 
basic helicity structure of the process under consideration. The general tensor structure 
of the amplitude is described in Section 3 and expressed through helicity amplitudes in 
Section 4. The calculation of the two-loop amplitudes, their renormalisation and infrared 
properties and their simplification are described in Section 5. The two- loop helicity ampli- 
tudes are obtained in a closed analytic form. We performed several non-trivial checks on 
the results, which are described in Section 6. We conclude with an outlook in Section 7. We 
enclose appendices with the analytical form for the one-loop and two-loop helicity ampli- 
tudes in the decay kinematics V — t- ggg and V — t- gg'j. The helicity amplitudes continued 
to the regions relevant for vector-boson-plus-jet and vector-boson-plus-photon production 
at LHC are given in Mathematica format together with the arXiv submission of this paper. 

2. Kinematics and notations 

The production of a massive vector boson V = (Z^,^*) and a gluon (photon) through 
gluon-gluon fusion is related by crossing to the decay of a massive vector boson to three 
gluons (two gluons and a photon) and has the same kinematics as vector-boson-plus-jet 
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production qq — t- Vg, qg — t- Vq and vector-boson-plus-photon production qq — t- Vj. Tech- 
nically the calculation of the two-loop QCD corrections to the gg — > Vg and gg — >■ V'j 
amplitudes is thus similar to previous calculations for 3j-production, vector-boson-plus- 
photon production and H 3 partons, which have been derived to two-loop accuracy in 
QCD [18-20]. 

In the following we will focus on the decay kinematics, while the crossings relevant for 
1^-plus-jet and IZ-plus-photon production at hadron colliders will be discussed in section 4.3. 
The relevant partonic subprocesses are: 

I'iPb) + 1^{V&) V{q) g{pi) + g(p2) + gips) , 

r{p5) + l^ipe) ^ V{q) ^ g{pi) + g{p2) + 7^) , (2.1) 

where we included the production of the vector boson V through lepton-antilepton anni- 
hilation. 

In the framework of massless QCD interchanging the virtual photon with a Z boson 
amounts only to a proper re-weighting of the final result. Moreover, note that we always 
assume massless fermions in the initial or final state. 

The momentum of the vector boson is given by 

g/^=pA*+p^+p^. (2.2) 
It is convenient to define the usual invariants 

which fulfil 

(f' = {P\ +P2+ Pzf = Sl2 + Sl3 + S2Z = S123 , (2.4) 

as well as the dimensionless invariants 

a; = si2/si23, y = si3/si23, = S23/S123 , (2.5) 

which satisfy x + y + z = 1. 

In the decay kinematics V — >• ggg/gg'J, as in the 3j case, q^ is time-like (hence positive) 
and all the Sij are also positive, which implies that x,y,z all lie in the interval [0; 1], with 
the above constraint x + y + z = 1. 

The helicity amplitudes can be expressed as a product of a partonic current and a 
leptonic current L^: 

MP5,P6;gi,g2,b3) = L>'{p^\pQ)S^{gr,g2-M) (2-6) 

where gi = g{pi), and 63 = 6(^3) labels a generic massless gauge boson. In our case b = g,'-f 
in V ggg and V — )■ ggj respectively. 

The purely vectorial tree-level leptonic current reads: 

L^{pb,P&) = v{p6)j''u{p5), (2.7) 
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where in the case of an incoming lepton-antilcpton pair L^j^^p^ jPe ) corresponds to a left- 
handed current, and Li^{p^ ,Pq ) to a right-handed current: 

Ll{p^,Pt) = v^{p^)ru.(jP^), L''n{pt,Pe) = v^{pe)l^ n+{p,). (2.8) 

Only the partonic currents receive contributions from QCD radiative corrections, and they 
can be perturbatively decomposed as: 

5^(9i;52;93) = V^sd-''''''' [(t) 4'^ (51; 52; 53) + {^fsj^\gr,g2;93) + 0{al) 

S^{9i;g2;i3)= V4^<5«i<^^ [ (^) 4'^(5i;52;73) + (^)' -^f (5i;52;73) + 

where we factored out the overall colour factors S""^""^, (/"^"^as 
The general form of the gauge boson coupling to fermions is: 



yV,/i/2 ^ _^gr^,/i/2 ^ith e = V47ra, (2.9) 

whose explicit form depends on the gauge boson, on the type of fermions, and on their 
lidicitjiGS * 

= 4./. {^) + 4./. 7„ (ii^) . (2.10) 

The left- and right-handed couplings are identical for a pure vector interaction, and are 
in general different if vector and axial- vector interactions contribute. Their values for a 
photon are 

^}lf. = <f. = -'h^hh^ (2-11) 

while for a Z boson 



z it - sin^ ^wef, ^ z single/, 

^f^f^ - sin0,„cos0,„ ^f^f' ' ^f^f^ - ^f^f' ■ ^^-^^^ 



The vector boson propagator can be written as: 

where A|fj^(g,^) and Dv{q) are, respectively, the numerator and the denominator in the 
gauge: 

KA<1,0 = {-9,. + (1 - 0^^^) , (2.14) 

Dz{q) = [q^ -Ml + iVzMz), (2.15) 

Dril) = (2-16) 

In the narrow- width approximation we can simplify expression (2.15) to 

Dz{q)^iTzMz and q^ = (2.17) 
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where Mz is the mass of the Z boson, while Tz is its decay width. 

Since we do not consider any electroweak corrections, the vector boson V is always 
coupled to a fermion line which allows us to neglect the dependence (or equivalently 
to put ^ = 1). A further consequence is that the total amplitude is proportional to the 
charge weighted sum over the quark flavours, such that all electroweak couplings can be 
collected into a multiplicative factor Qy. With this notation we obtain for an incoming 
right-handed lepton-antilepton pair, for the different choices oiV = (7*,^), and helicity 
configurations (/ii, /i2, /is): 

■My(P5+,P6 ;gl^g2^^^) = -^(4vr«) ^^5"^ , 4'^'^'^^(P5,P6;ffi,g2,b3), (2.18) 



Dv{pb+P&) 



In case of F = 7* we find 



Q^*=Egeg, (2.19) 
Ql,*=T.qel (2.20) 

where the sum runs over the quark fiavours in the loop. 

In the case oi V = Z we have a contribution from the vector component of the Z 
boson, which is given by 

Q| = ^Eg(^^. + 0, (2.21) 
Ql = lE,(^g^g + ^g>., (2.22) 

but also a contribution involving its axial coupling. This contribution vanishes for Z — )■ gg'^ 
due to charge conjugation invariance, already before summing over the quark flavours in 
the loop. On the other hand, in the case oi Z ^ ggg it vanishes only after summing over 
the quark flavours. 



3. The general tensor structure 

In order to extract the helicity amplitudes from a generic QCD process different approaches 
can be attempted. One possibility is to decompose the amplitude into linearly independent 
tensor structures, whose number and form are entirely determined by symmetry consid- 
erations and which are completely independent on the loop order we are interested in. 
The entire loop-dcpcndcncc is then contained in the scalar coefficients which multiply the 
relevant tensor structures. In order to single out these coefficients we apply projectors 
defined in d-continuous dimensions directly on the Feynman-diagrammatic expression for 
the amplitude [18-21]. 

Using Lorentz invariance one can show that there are 138 independent Lorentz struc- 
tures which can contribute to the partonic current [21]: 

S^'''P'' = aig^^'g^^ + a2g^'^g'"' + asg^^g^^ 
iij2=i 
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3 

+ (^hhJiuPj.PhPjsPjr (3-1) 

jl,J2J3,J4=l 

Not all these tensors will be relevant for our computations. Defining the physical amplitude 
contracted with the external polatization vectors of the three massless on-shell bosons: 

SiM{gi\g2-M) = Sij,„pa{pi;p2;P3) ei{g) e^{g) e'^{b), (3.2) 

we see that many of the structures do not contribute because of the transversality condition: 

Ci ■ Pi = 0, with i = 1,2,3. 

This reduces the number of independent tensors to 57. One way of proceeding is then to 
apply Ward identities for the massless bosons 

gl^-ypa ^ gt^upa ^p ^ gf^upa ^ q ^3 3^ 

which lowers the number of relevant structures down to 18. Applying finally current con- 
servation for the massive boson 

S'^'P'' e\ 4 el p^ = 0. (3.4) 

further reduces the number of independent tensor coefficients to 14. 

By requiring the amplitude to be invariant under the exchange of the three (two) gluons 
one can find further relations among these 14 coefficients with interchanged arguments. 
This allows to perform different checks on the final result (see section 6). 

Once the tensor structure is known, one can compute d-dimensional projection opera- 
tors that applied on Sfxvpa extract each of the 14 coefficients. The tensors and the projectors 
contain a large number of individual terms. Therefore applying them to an amplitude in a 
Feynman-diagrammatic approach will in general result in a large number of contractions 
with a huge proliferation of terms. 

Moreover, it must be noted that the basis of tensors is not unique, namely that any 
set of 14 tensors, obtained as independent linear combinations of those found above, can 
be chosen. Choosing suitable linear combinations of the above tensors can simplify their 
structure substantially. 

For all these reasons we decided to follow a simplified approch, which nevertheless al- 
lows us to retain the full information on the process. It is well known that when performing 
a computation with a large number of external bosons a specific gauge choice can highly 
simplify the intermediate steps of the calculation, while gauge invariance ensures that the 
final result for the amplitude must be independent on the choice made. Following this idea, 
instead of imposing gauge invariance on the tensor structures, we chose to fix the gauge of 
the external particles in order to symplify the tensor structures as much as possible. 
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Naively one would expect the loss of gauge invariance on the tensors, together with 
the loss of part of the symmetry due to the gauge choice performed, to be a drawback 
of this approach. However, one can show that once these 14 coefficients are known, the 
full gauge-invariant tensor can be reconstructed. In particular one can find linear relations 
among the 14 coefficients obtained imposing the gauge fixing and the 14 coefficients of the 
gauge invariant tensor, as outlined in the following section. 

3.1 The gauge- fixed tensor structure 

Following the above reasoning, we replace the condition (3.3) with a gauge choice on the 
external on-shell bosons: 

ei • P2 = €2 • P3 = €3 • Pi = 0. (3.5) 

This choice is arbitrary and could be substituted by any other set of gauge conditions. The 
advantage of this particular choice is to produce extremely compact tensor structures. 

Fixing the gauge of the external bosons reduces the number of independent tensors to 
18. Also in this case we impose current conservation (3.4) on the and end up again with 
14 tensor structures. As expected, the number of independent tensor structures obtained 
in this way is the same as for the gauge-independent tensor. 

We decompose the parton current as 



14 



(3.6) 



where the coefficients are functions of the mandelstam variables Af^ 
and their explicit values differ in general if 6 is a gluon or a photon. 
Finally, the gauge-fixed tensors read: 



Af'\si2,Si3,S23) 



= ei • Pa ea • p2 - 62 ■ pi ea ■ P2 el, 



ei • Pa e2 ■ Pi - £2 • Pi ea • P2 e^, (3.7) 



= 61-62 

?t = ei • 62 
= ei • €2 



ea • P2 Pi 



ea •P2P2 



ea •P2P3 



(Sl2 + S13) 



{Sl2 + S23) 



(sia + S23) 



(3.8) 
(3.9) 
(3.10) 



ri' = ei • ea 



T!^ = ei • ea 
= ei • ea 



e2 • Pi Pi 



e2 • Pi P2 



e2 -piPa 



(Sl2 + '^13) 



(Si2 + S23) 



(sia + S23) 



(3.11) 
(3.12) 
(3.13) 
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Tg" = €2 • 63 
TiQ = £2 • £3 



(Sl2 + S13) 



ei -^3^2 



n jsu + ^23) 



(gl3 + ^23) fi 



(3.14) 
(3.15) 
(3.16) 



T12 = 62 • Pl 63 • P2 
Ti'3 = 62 • Pi 63 • P2 
^14 = €2 • Pi 63 • P2 
l(b) 



yit (gl2 + 513) M 

ei -psPi 2 ^1 



ei •P3P2 



ei -PsPs 



(gi2 + ^23) ij. 

_ (Sl3 + ^23) fi 



(3.17) 
(3.18) 
(3.19) 



The relations among the ' and the coefficients of the gauge invariant tensor can be 
found by performing on the latter the gauge fixing (3.5). This procedure obviously does 
not affect the scalar coefficients which multiply the tensor structures. One ends up then 
with 14 new tensor structures which can be related through hnear combinations to those 
obtained fixing the gauge from the beginning. In this way the gauge invariant tensor can 
be fully reconstructed. We have verified this procedure by comparing our one-loop result 
with the literature [4] where the results are given for an on-shell Z boson, and a different 
gauge choice is used (see section 6). 

Once the tensor structure is known, one can obtain the coefficients Af^ by applying 
a set of projectors Pn{Af'^) on the Feynman-diagrammatic expression of the amplitude 
defined such that 



Y,P'{Af^)S^{pi,P2,P3) = AfK 



spin 



Note that the projection has to be performed in d dimensions, and that special care 
has to be taken in performing the polarization sums when applying the projectors on the 
single diagrams. In particular one has to consistently use a physical polarization sum which 
respects the gauge choice (3.5): 



spin 
spin 

E^rb3)e^(P3) 



Pi •P2 

...p^Ps+pM 

P2 •P3 



spin 



P3 -Pl 



(3.20) 
(3.21) 
(3.22) 
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The projectors themselves can be decomposed in the tensor basis and take the form: 

14 

P^iAf'') = XiiAf^) T*^ (3.23) 

where ih.e Xi{Af) are functions of d and the kinematical invariants Sjj. 
4. Helicity amplitudes 

By fixing the hehcities of the external massless bosons the partonic current can be cast in 
the usual spinor helicity notation [22]. There are two independent helicity configurations 
in the gggV-case, and three independent helicity configurations in the ggjV-case, from 
which all the others can be obtained. In the following we discuss separately the two cases. 

4.1 gggV: The amplitude in spinor helicity notation 

We start off considering the gggV-case. We choose as two independent helicity configu- 
rations {gt , g2 : 9^) {di J dt ^ dt)- order to include the spin-correlations with the 
leptonic decay products we contract the partonic current with the leptonic current for 
fixed helicities of the initial state Icptons. This also helps to further simplify the result. 
Consider the production of the vector boson V through lepton-antilepton annihilation: 

The leptonic currents (2.7) are 

L'kiPt^P^) = [6 It'^I 5), L^LiP^,pt) = [5 It'^I 6) = [L'^M^Pe)]*- (4-1) 

Performing the contraction and making use of Schouten identities and momentum 
conservation we end up with: 

A^r''\p5,P6; 91, 92,93) = L''^{pt,PQ)S^{g+,g2,g2) = -j= ^^2)(13)[2 3] 



X <^ (2 5) (3 5) [5 6] ai (x, y, z) + (2 3) (2 5) [2 6] ^2 {x, y, z) + (2 3) (3 5) [3 G]a3{x,y, 



z 



(4.2) 



A^R^^\p5,P6;9i,92,93) = L'^{pt,PG) Si^igf, 92^93 ) = 



I [13](15)[16] [23]{25)[26] [23](25)[16] 1 

\ (12)(23) ^^^'2/'^^+ (12)(13) /^2ia;,2/,zj+ ^^2){23) 

(4.3) 

where the coefficients aj and Pi are linear combinations of the 14 tensor coefficients A^. As 
an explicit example we write down the relations for the aj: 



ai{x,y,z) = -{si2 + Sis) ^2 + ^9 + ^^12 - (si2 + S23) 



Al + ^10 + ^^13 
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- {Sl3 + S23) 

a2ix,y,z) = -S12 

- (Sl3 + ■523 
a3{x,y,z) = Si3 



An + fAu 
A2 + Ag + ^Ai2 

An + f Ai4 



(Sl2 + Si3 + S23) 



^1+^10 + ^^13 



A2 + Ag- An + ^^12 - '-fAu 



(4.4) 

(4.5) 
(4.6) 



The corresponding relations for the are slightly longer and we do not reproduce them 
here for brevity. There are in total 16 different helicity configurations. Prom the above 
expressions for \p5,P6; 9i, 92, 93) and A^^'^'^\p5,pe; gi, 92, 93), all the other helicity 

amplitudes can be obtained by parity conjugation and permutations of the external legs. 
We find: 



L'niPt 


Pq)Sh{9i 


,92 


93 ) 


= 4- 


'^~^{pb,P&;9i,92,93) 


= 4^- 


~\p5,Pe;92,9i,93) , 


LUpt 


P&)Si,{gi 


,92 


9t) 


= 4- 


~^\P5,P6;91,92,93) 


= 4- 


~\p5,P6;93,92,9l) , 


L'M 


P6)Sf^i9t 


,9t 


93) 


= 4* 


'^\pb,P&;9i,92,93) 


= [At 


\p6,P5; 93,92,91)]* , 


L'M 


P6)Sf,igt 


92 


9t) 




~^\pb,p&; 91,92,93) 


= [At 


~\P6,P5;92,91,93)]* , 


LUpt 


P6)Sn{9i 


,9t 


9t) 




^^\P5,P6;91,92,93) 


= [A't 


'\P6,P5;91,92,93)]* , 




P%)Sf,{9i 


,92 


93) 


=4- 


~~\pb,p&; 91,92,93) 


= [4+++)(p6,P5;5i,52,53)]*- (4.7) 



The corresponding amplitudes for right-handed leptonic current can be obtained by simply 
interchanging ■<-> pe- Note that the complex conjugation operation has to be applied 
only on the spinor structures in (4.2) (4.3), and not on the coefficients aj, Pj. 

The unrenormalised helicity amplitude coefficients are vectors in colour space and have 
perturbative expansions: 



as 
2^ 



^^?^'"'^+(g)'f^r'"° + o(«^) 



(4.; 



for Qg = ai,Pi. The dependence on {x,y,z) is again implicit. 
4.2 gg^V: The amplitude in spinor helicity notation 

In the gg-jV-case there are three independent helicity configurations. Two of them can be 
chosen identical to those in the gggV-case, namely {gf , 92 ,'^3) and {gf , g2 ,13), the third 
is taken as (5^ , ^2" > 73" ) • 

Fixing the helicities and contracting with the right-handed lepton current we have: 



A 



R \p5,P6;gi,g2,-f3) = L'^{pt,P(i)SM,92 ,13 ) = ^^2)(i 



(2 3) 



3) [2 3] 



X <^ (2 5) (3 5) [5 6] 771 {x, y, z) + (2 3) (2 5) [2 6] m (x, y, z) + (2 3) (3 5) [3 6] % {x, y, z) 



(4.9) 
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(12) (2 3) 



(12)(13) 



(12)(23) 



(4.10) 



[12] 



^(12) [13] [2 3] 

X J [1 2] [1 6] (1 5) n {x, y, z) + [1 2] [2 6] (2 5) T2 {x, y, z) + [1 6] [2 6] (6 5) (x, y, z) 



(4.11) 



From A^^ Hp5,P6; 92, gi,l3), ^^^^ Hp5,P6; 9i, 92,^3) and ^^^^''(^5,^6; 51, 52, Ta) 
all the other helicity configurations can be obtained by parity and charge conjugation: 

AR^~\p5,P6;gi,g2,73) = -4^ \p5,P6;g2,gi,j3) 
~'^Hp5,P6;9i,52,73) = [A^R'^~\p6,P5;gi,g2,73)]* 

A^R~^\p5,P6; 91, 92,13) = [a''r'"'\p6,P5; 92, 91,13)]* 
Ar^^\p5,P6; 91, 92,13) = [a'-^~~\p6,P5; 91, 92,13)]* 
AR^~~\p5,P6;9i,g2,i3) = [A^R^^\p6,P5;gi,g2,i3)]*- (4.12) 

As before, the left-handed helicity amplitudes can be found by the exchange p5 pe, and 
the complex conj ligation has to be performed only on the spinor structures and not on the 

coefficients r]j, 9j, tj. 

The unrenormalised helicity amplitude coefficients are vectors in colour space and have 
perturbative expansions: 



or = V4^(5«i"2 



|)l^«-^+(^)^i.?)- + (9(a3) 



(4.13) 



T • \2TrJ ■ V27r 

for rt-y = r)i, Oi, Ti- The dependence on {x, y, z) is again implicit. 

4.3 Analytic continuation to the scattering kinematics 

In order to compute the two- loop contributions to F-plus-jet and F-plus-photon produc- 
tion at hadron colliders, the helicity amplitudes must be continued to the appropriate 
kinematical situations. 

The relevant partonic subprocesses are: 

giPi) + 9{P2) ^ 9{-P3) + y{q) 9{-P3) + 1^{P^) + I'iPe) , (4.14) 

g{P2) + g{P3) ^ g{-Pi) + vio) ^ g{-Pi) + ^+(^5) + r (pe) , (4.i5) 
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where the second crossing is required to fully account for all helicity combinations, and 



9{Pi) + 9{P2) li-ps) + V{q) 7(-p3) + I'^iPb) + I (Pe) • 



(4.16) 



With the notation above the definitions of the helicity amplitudes in terms of momen- 
tum spinors (4.2) (4.3) and (4.11) remain unchanged under crossing. Considering in fact 
an outgoing leptonic current defined as: 



with 



we find that: 



L^{P5,P6) = u{P6) 1^ v{P5), 
out 



L'nipt^Pe) . =[6|7l5)=L^(p5-,p+ 
Ll(p^,Pt) . =[5 17^^16) = L^(p+,P6-; 



(4.17) 
(4.18) 



out 



out 



This means that the expressions for the helicity amplitudes defined in the two sections 
above remain unchanged provided that is now considered as the label of the antilepton 
and pe the one of the lepton. 

Special care has to be taken in the analytic continuation of the helicity coefficients 0,g 
and $7^. In the kinematical situation in (4.14) and (4.16) remains time-like, but only 
si2 becomes positive: 



or, equivalently. 



> , Sl2 > , Si3 < 0, S23 < 



X > , y <0 , z <0 . 



(4.19) 



(4.20) 



As shown in [23] (where this region is denoted as (2a)+) and used for example in [19], this 
kinematical situation can be expressed by introducing new dimensionless variables 



Ul 



Sl3 
Sl2 



y 

X 



Vl 



Sl2 



1 

X 



(4.21) 



which fulfil 



< «i < , < < 1 . 



To account for all helicity combinations in the case of gg — )• gV ^ also the kinematical 
situation (4.15) must be considered. In this case we have 



g > , Sl2 < , Si3 < 0, S23 > , 



(4.22) 



This can be treated with the following choice of variables [23] (this region is denoted as 
(4a)+) : 

si3 y I 



U2 



S23 



V2 



S23 



(4.23) 
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which fulfil again 

0<U2<V2, 0<V2<1 . 

Note that the two kinematical regions (4.14) and (4.15) are turned each other by the 
permutation pi 0^3, in particular one has: 

Ul{Pl ^ Ps) = U2 
Vi{pi O ps) = V2. 

As shown in (4.7), in the gggV-casc, in order to obtain all the different helicity configura- 
tions, we also need to exploit the Bose symmetry of the external gluons. It is now clear 
that whenever the permutation pi -f-)- p^ is performed, this only amounts to switching from 
region (4.14) to region (4.15). 

We provide the one-loop and two-loop coefficients in all relevant regions in Mathemat- 
ica format together with the arXiv-submission of this paper. 

5. Outline of the calculation 

The two-loop corrections to the coefRcients can be evaluated through a calculation of 
the relevant Feynman diagrams. The calculation proceeds as follows. The diagrams con- 
tributing to the process are produced using QGRAF [24]. In the gggV-case there are 12 
diagrams at one loop and 264 at two loops, while in the ggjV-case there are 8 diagrams 
at one and 138 at two loops. The tensor coefficients are evaluated analytically diagram by 
diagram applying the projectors defined above. As a result, one obtains the tensor coeffi- 
cients in terms of thousands of planar and non-planar two-loop scalar integrals, which can 
be classified in two auxiliary topologies, one planar and the other non-planar [25]. In order 
to do so, one needs to perform both shifts in the integration variables and permutations on 
the external legs. All the routines needed for this purpose have been coded in FORM [26] 
and checked against the new automated shift-finder implemented in Rcduze2 [27] . Through 
the usual IBP identities [28] one can reduce independently all the integrals belonging to 
these two auxiliary topologies to a small set of master integrals. This reduction is per- 
formed using the Laporta algorithm [29] implemented in the Reduze code [27, 30] . All the 
masters for the topologies above are known as series in the parameter e = (4— d)/2 through 
a systematic approach based on the differential equation method [31,32]. The masters are 
expressed as Laurent expansion in e, with coefficients containing harmonic polylogarithms 
(HPLs, [33]) and two-dimensional harmonic polylogarithms (2dHPLs, [31]). Numerical 
implementations of these functions are available [34]. For all the intermediate algebraic 
manipulations we have made extensive use of FORM [26] and Mathematica [35]. The 
two-loop unrenormalised helicity coefficients fi^^^'"" can then be evaluated as linear com- 
bination of the tensor coefficients. The whole computation is performed in the euclidean 
non-physical region, where the amplitude is real. The final result is then analytically con- 
tinued to the physical regions relevant for Z + jet/j production at LHC, as thoroughly 
discussed in [23] and in section 4.3. 
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5.1 UV Renormalisation and IR subtraction 



Renormalisation of ultraviolet divergences is performed in the MS scheme by replacing the 
bare coupling uq with the renormalised coupling = as{n'^), evaluated at the renormal- 
isation scale iJp. Since there is no tree level contribution to the amplitude, we only need 
the one loop relation between the bare and renormalised couplings: 



where 



(5.1) 



= (47r)'e"^^ with Euler constant 7 = 0.5772 . . . 

and /iQ is the mass parameter introduced in dimensional rcgularisation to maintain a di- 
mensionless coupling in the bare QCD Lagrangian density. /3o is the first coefficient of the 
QCD /3-function: 

^^^ llC.-4T,iV. _ (5.2) 

with the QCD colour factors 

N"^ - I 1 
Ca = N, Cf = ^^, Tr = -. (5.3) 

The renormalisation is performed at fixed scale 11^ = q^. The renormalised helicity coeffi- 
cients read: 



0« = ^^if^W'"", 



After performing ultraviolet renormalisation, the amplitudes still contain singularities, 
which are of infrared origin and will be analytically cancelled by those occurring in radiative 
processes of the same order. Catani [36] has shown how to organise the infrared pole 
structure of the one- and two-loop contributions renormalised in the MS-scheme in terms 

of the tree and renormalised one- loop amplitudes. The same procedure applies to the tensor 
coefficients. Since there is no tree level process contributing, their pole structure can be 
separated off as follows: 

^^f^ = 4'^(^)^^^'^+^^f^''°"^ (5.6) 

where again 6 = g', 7. 
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In the two cases the operator I^^ (e) is given by 



4"(. 



=£7 



2r(l - e) 



1 . A 



(Si2 + Si3 + S23) 



1 



A^e 



312 



(5.7) 
(5.8) 



where, since we have set //^ = S123: 



Sij — 



S123 



(5.9) 



Note that on expanding Sy, imaginary parts are generated, depending on which kinematical 
configuration we are working in. In the decay kinematics Z — t- ggg / gg^ for example we 
have that ah the Sij become positive, so that all three terms will generate imaginary parts 
whose sign is fixed by the small imaginary part +iO of Sij. On the other hand if we are 
interested in the scattering kinematics gg — > Zg / Zj only S12 or S13 become positive, with 
the usual Sij + iO prescription. 

For the infrared factorisation of the two- loop results, the renormalised one- loop he- 
licity amplitude coefficients are needed through to 0{e'^). Their decomposition in colour 
structures is straightforward, namely the whole colour dependence is in the overall factors 
^010203 g^j^^ ^0102 fQj. gggV aud ggjV respectively. 



an^{x,y,z) 



(5.10) 



The expansion of the coefficients through to yields HPLs and 2dHPLs up to weight 4. 
The explicit expressions are of considerable size, such that we only quote the e^-terms in 
the appendix. To this order, the coefficients had been derived previously [4] in terms of 
logarithms and dilogarithms. The expressions through to 0{e^) in Mathematica format 
are appended to the arXiv submission of this article. 

The finite two-loop remainder is obtained by subtracting the predicted infrared struc- 
ture (expanded through to 0(e°)) from the renormalised helicity coefficient. We further 
decompose it according to the colour structures as follows: 



NAn, + ^Bn, + NfCn,., 



(5.11) 



The helicity coefficients contain HPLs and 2dIIPLs up to weight 4. The size of each 
helicity coefficient is comparable to the size of the helicity-averaged tree times two-loop 
matrix element for 3j production quoted in [25] , and we decided not to include them here 
explicitly. The complete set of coefficients in Mathematica format is attached to the arXiv 
submission of this article. 



5.2 Simplification using the Symbol formalism 

After the computation of the amplitudes and subtraction of UV- and IR-divergences we 
used an in-house implementation of the algorithm described in [37] to express the result as 
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far as possible in logarithms and poly logarithms of functions of the kinematic invariants. 
The GiNaC libary was used to evaluate the 2dHPLs [38] and the implementation of the 
PSLQ algorithm contained in the arprec library [39] to find the parts mapped to zero by 
the symbol map. 

It is well known that up to transcendental weight three all two-dimensional harmonic 

polylogarithms can be expressed this way. However, for weight four this is not always 
the case. In [40] it was conjectured that a combination of 2dIIPLs can be expressed in 
logarithms and polylogarithms if and only if its symbol fulfills a certain symmetry condition. 
In the present case, we found this condition in general not to be fulfilled and were also not 
able to express our result in logarithms and polylogarithms only. Nevertheless we reduced 
the number of required functions in all kinematic regions as far as possible, having to resort 
to 17 2dHPLs of weight four. 

In the past, surprising relations between certain QCD and N=4 SYM amplitudes 
have been found, for example in the case of if — )• ggg at two loops in the heavy-top- 
limit [20,41,42]. In the leading color part of the finite two- loop amplitude, the weight four 
contribution without a rational factor was found [41] to be helicity-independent and equal 
to the three-point form factor remainder function in planar N=4 SYM. In the present cases, 
however, no such relation could be observed. This feature can be understood from the fact 
that, in contrast to the Higgs amplitudes, no purely gluonic contribution is present here, 
due to the internal quark loop coupling to the vector boson. 

6. Checks on the result 

Several non-trivial checks were applied to validate our results. 

1. As a first check we computed all 14 tensor coefficients in (3.6) at one-loop order for the 
gggV-case, and we verified that we can reproduce the results in [4] up to order O(e^). 
Performing this check was not entirely trivial. In [4] the results for the one-loop 
helicity amplitudes are given in the case of an on-shell Z with a fixed polarization. 
Moreover, the amplitudes for different helicity configurations are given choosing an 
explicit representation for the polarization vectors of the external particles. This 
representation docs not respect the gauge choice performed in (3.5), so that we cannot 
naively start from our tensor structure and fix the polarization vectors in the same 
way to reproduce their result. Nevertheless, as explained in section 3, the full gauge- 
invariant tensor can be fully reconstructed taking suitable linear combinations of 
the tensor coefficients of the gauge-fixed tensor. Once the gauge-invariant tensor is 
known, one can then use the explicit representation of the polarization vectors given 
in [4] and demonstrate the analytic agreement of the expressions. 

2. We computed all the 14 tensor coefficients both at one-loop and at two-loop order, 
in the gggV- and in the gg^V-case. Following the procedure outlined in section 3, 

we obtained the 14 coefficients of the gauge invariant tensor for both processes, and 
we verified that they respect the expected symmetry relations under permutation of 
the external gluons. 
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3. The IR singularity structure of our results agrees with the prediction of Catani for- 
mula [36], see section 5.1. 

4. We compared the helicity amplitudes 0.^^ for the gggV- and the gg^V-ca^se. We 
verified the following identities for the one-loop amplitude coefficients: 

2aaj{x,y,z) = ari^{x,y,z), 

'^a/3jix,y,z) = a0.{x,y,z). j = 1,2,3. (6.1) 

5. Finally, we performed the same comparison at two-loop order, finding: 

2Baj{x,y,z) = Br,^{x,y,z), 

2Bp^{x,y,z) = Bg^{x,y,z), i = 1,2,3, (6.2) 

which follow from the structure of the underlying two-loop diagrams. The subleading 
colour coefficients B are unaffected by renormalisation and infrared subtraction. No 
relation of this type can be found for the coefficients Cn^, which are determined 
purely from renormalisation counterterms and IR subtraction, which differ in the 
cases 6 = 5, 7. 

7. Conclusions and Outlook 

In this paper we presented the two-loop corrections to the helicity amplitudes for the pro- 
cesses gg Vg and gg V^. We performed the calculation in dimensional regularisation 
by applying d-dimensional projection operators to the most general tensor structure of the 
amplitude. We showed how an explicit gauge choice can reduce considerably the com- 
plexity of the basic tensor structures appearing while retaining the full information on the 
gauge-invariant amplitTides. We expressed our results in terms of dimensionlcss helicity 
coefficients, which multiply four-dimensional spinor structures. Wc extracted the infrared 
singularities by means of an infrared factorisation formula and provide compact analytic 
expressions for the finite part of the two-loop helicity coefficients in all relevant kinematical 
regions. 

The matrix elements derived here contribute to the NLO corrections to the gluon- 
indiiced production of and Z + j final states at the LHC. Viewed in an expansion in 
the strong coupling constant, these contributions arc formally N^LO as far as the reactions 
pp — F7 + X, pp ^ Vj + X are concerned. However, due to the large gluon-gluon 
luminosity at the LHC, these contributions could be comparable in size with the NNLO 
corrections to qq — > Vg, qg — > Vq and qq — > V^. Their inclusion will also help to stabilise 
the substantial scale dependence of the gluon-induced subprocesses, which were known only 
at Born-level up to now. 
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A. One-loop helicity amplitudes 
A.l V ggg at one loop 

We reproduce the leading order O(e^) for V ggg. The complete expressions up to order 
O(e^) can be found in the attachments of the arXiv version of this paper. 



oai ix,y,z) 



2x I ) log(x) + 2y (-^ ) log(y) 

1 — X z J \\ — y z ' 



{I - x)x + {I - y)y 



— + log(x) log{y) 

- ( log(l - x) log(x) + Li2(x) ) 

- ( log(l - y) log(y) + Li2 (y) ) 



(A.l) 



aa2{x,y,z) = 2y 



1 1\ 2x{2y + z) 2xy{z + {1 - y)i2y + z)) ^ 



1-y z 
x (2y2 + 2yz + z^) 



(l-y)2^2 



vr 



-I- log(a;) log(y) 



- (log(l - x) \og{x) + Li2(x) ) 

- ( log(l - y) log(?/) + Li2(y) ) 



(A.2) 



aaz{x,y,z) = -aa2{y,x,z) , 

afii{x,y,z) = -2^1-^^ , ai3^{x,y,z) = -2(^1-^^ , 
a/^six, y,z) = - 4. 

A.2 V — > gg^ at one loop 

At one loop, the two amplitudes are related to each other as follows: 

ar,^{x,y,z) = 2aa^{x,y,z) fori = 1,2,3, 
{x,y,z) = 2afj^{x,y,z) for z = 1, 2, 3 , 

On {x, y,z) = 2 {z, y,x), {x, y,z) = 2 a^^ {z, y, x) , 
aTz{x,y,z) = 2aai{z,y,x) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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B. Two- loop amplitudes: all-plus helicity coefficients 

Due to the length of the resulting expressions we chose to reproduce only the all-plus 

{Qi 1 dt ■• dt / ^t) helicity amplitudes of both processes in the decay region, which are con- 
siderably shorter than the other helicity combinations and contain only functions up to 
transcendental weight two. The full result can be found in the attachments to the arXiv 
submission of this paper in Mathematica format. 

B.l F — > ggg at two loops 

The coefficients for the {gt , gt ^ 9t) helicity configuration are: 



^ f r,f 1 1 1 \ 4z 4^2 



X x^ 



1 (2,{l-y) 2(-l + 2y) , 3(1 -y) , U{l-y)z Hz^X ^ 

H 1 5 5- TT 



12 \ xy y^ yz y^ y^ 

If 6 6 A2x + lly\ , , , 1 / 47\ , , , 

1 -42 + 42x + 31y 6 6 \ , , , 

- J ( - + — ) log (x) log (y) - J ( - + — ) log(y)log(z) 

2 \y yzj 2 \y xy J 

7(1 - x)x 1 7x 1\ , , , , , , 

y6 yl yl yj 

^ 1 /2(1 - 7z) _ 3 ^ 3(1 -z) ^ 14(1 - z)z 



2 V y2 z yz y^ 

X (log(l - x) log(x) -I- Li2(a;)) 

+ + — + — ) (log(l - 2/) log(?/) + Li2(y)) 

1 fUx'^ 142;(l-y) 3(1 - y) _^ -2 + 3y 



2 \ y^ y^ xy 

X (log(l-z)log(z)+Li2(z)) , (B.l) 



Ap2ix,y,z) = - — 



1 

27 \ 1 

4x 



27 



X 



+ 



27 Ay Ay 



+ 
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y 



O "1" ^ I o Q 



X 



X^ 



Ay 



y^ 



1 



1 ] ITT 

Z. 



_ 54\ _ 11 

y \ z-^ ) z ) 2 

1 /rr2(14 + y(-58 + 45y)) ^ lAx^{-2 + 3y) 2(1 - y)y(-8 + 21y) 



12 \^ y^z^ y^z^ xz^ 

4x(l-y)(-4 + 5y) , 3(2 + 3y(-4 + 5y)) , 14(1 - y)V , 14x 



yz^ 



+ 



+ 



2 2 



+ 
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1 f53x 42x2 (_47_53(_2 + x)x)?/\ , 
6 \ z yz (1 — xj^'z J 

^1 f 53y 42y2 ^{47 + 53i-2 + y)y)\ 

+ 6 [-— o^^F^ J ^^^^^ 

1 / 42x 42y 11\ , , , 

-- 31 + + ^ + — log(z) 

D \ y X z J 



1 /x2 x(2-12y) y(2 + y) 

2 (^^2 + -2 + ^2 ) log(^) log(y) 

14(1 _ 2(l-z)(-l + 7z) 

1 / 14y2 -2 + 14y + ^ 

2 ' + ^ + 



o 1 1 + — n 7z ) log(^) log(^) 



l + ^+ "^"?^"- llog(y)log(.) 



^ 1 / 14(1- + IV _ - ^) _ 2(l--;)(-l + 7z) 

2 \ y2 

^ 3 + 2(-2 + z)z ^ (log(l - X) log(x) + Li2(x)) 
1 / 14(1 - zf 14x2 14x(l - z) 2(1 - z)(-l + 7z) 



2 V iC^ XZ 



^ 3 + 2i2 + z)z J ^^^^^^ _ ^^^^^^ ^ ^.^^^^^ 



z 

1-x 7x'^ ^ x(-8 + 7x) 7y 7y2 1 



(1 — x)y X x^ (1 — x)xz^ 

X (log(l-z)log(z) + Li2(z)) , 

(B.2) 



A^3(x,y,z) = --(^81-^-- + ^^^ + ^(^l-z + ^) I-IUtt 



1/5 14(1 - x)x -5 + 14x 14(1 - x)y 

TT: 3 ;5 h 



12 \ X y2 y X 

14y2 14(1 -x)x 5-14x\ ^ 

TT 



X^ 

1 / ^ 3 21(1 -x)x\ , , , 

3 V 1 - y a;^; / 

1 3 21(1 - z)z\ , , , 

1 / 14x2 14x(l - z) -h + z 



2\ z' z' )log(x)log(y) 

1/ U{l-z)z , -5 + 14z \, . X, . X 
~ 2 1^ ^ + y ' ^^^^^ ^^^"^^ 
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+ ^ (^2 - -^^2^ (14(1 - xfx - (1 - x)(-5 + 42x)y + (-5 + 42x)/)^ 

X (log(l - x) \og{x) + lA2{x)) 
+ \{^- ^^2^ (14(1 - yfy - (1 - y)(-5 + 42y)z + (-5 + 42y)z2)^ 

X (log(l - y) \og{y) + Li2(y)) 



1/14x2 14x(l - y) _ 14(1 - y)y -5 + 2y -S + HyX 

2^2/2 y2 3,2 y X ) 

X (log(l-z)log(z) + Li2(z)) , 

(B.3) 

N 1 1 1 - 

BpAx,y,z) =- 1 h 



1 — X y 1 — z 

I (2 l-y l-y 2(1 - y)^ , 2z^\ _2 



12 Vy y^ ^ 

+ ( - J) i°g(-) - - (T^) i°g(-) 

- \ log(x)log(y) - log(a;) log(z) 

(^)iog(y)iog(^) 

+ i ( — + ^ ) (log(l - ^) log(^) + Li2(a;)) 
+ ^ ( — + — ) (iog(i - y) log(y) + Li2(y)) 



2\yz xy 

U-^--] (log(l - ^) log(^) + Li2(^)) , (B.4) 
2 V y xy; 



1 1 2 

Bi3^{x,y,z) =- Tj h - 

1 — xl — 



1 , 2x2 , 2x , 2j/ , 2y^ , 1 - 2(1 - x)x 2(1 - x) , 

~T^ "^^ 2"^ ' ^^'^ 2 I ^ 

12 \ y X x^ 2; 

2x X x^z \ ^ ^ / X y \ , , , 
+ vT log(^) - 1 + - + - log(2) 



2 



(1 — x)2 z (1 — x)2y y y y X 

+ (i-(I^ + »(-i(T^ + i))M.) 

1 y2 \ , , , , 1 / 2x2 2x 
- o — + ^ log X log y - - 1 + _ + — 

2 2 \ y^ y 



1 / 2y 2y2 . ^ ^ , ^ ^ 
^ + X + ^^^^^ "^^''^ 
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/_^2x{l- z) _^ 1-2(1 -x)x 2(1 -X 



2 \ z 

X (log(l - x) log(x) + Li2(x)) 
_^ 1 2y _^ 2y2 ^ 1-2(1 2(1 - x 



2\ ' X ' x^ z^ z 

X (log(l - y) log(j/) + Li2(y)) 

+ ( 1 + ^^^^^^2^?-^) aog(l - ^) log(^) + Li2(^)) , (B.5) 



1/1 2{l-x)x , -l + 2x 2(1 
= -3- tt: 3-- -2 + 



12 \ x y x^ 

V 2(1 l-2x \ 2 

X X x\ (y xy , , . ^ 
log(x) + - + 7- ^ log(y) 



1 — X y z) \x {l — y)z 

(1 — X 1 z\ ^ , , 

+ + - log^ 

\ y 1-z xj 

1 ( (1 -x)x {l-y)y \ 



1 / 2{l -z)z , -l + 2z ^ M / N 

"2V y2 — + — - — ) l°g(^) l«g(^) 



+ ^ ( 2 - ^ ((1 - x) (2(1 - xfx + y + x(-7 + 6x)y + (-1 + Qx)y'^))\ 

^ \ y ^ J 

X (log(l - x) log(x) + Li2(x)) 
+ ^ (2 - -1^ ((1 - y) (2(1 - + z + yi-7 + 6y)z + (-1 + 6^)^^))^ 

X (log(l - y) log(y) + Li2(y)) 
^ 1 /2x2 _ 2x(l - y) _ 2(1 -y)y ^ -l + 2y ^ -1 + 2y \ 
2 \ x^ X V J 

X (log(l-z)log(z)+Li2(z)) , 

(B.6) 

C^i (x, y, =^1 - (3m - log(x) - log(y) - log(z)) , (B.7) 

C^,ix,y,z) =i (^1 - (3i7r - log(x) - log(y) - log(z)) , (B.8) 
2 

Cft (a;, y, 2) =2 (3i7r - log(x) - log(y) - log(z)) . (B.9) 
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B.2 V ggj at two loops 

The coefficients for the {gt ^92 ^it) heUcity configuration are as follows: 

2 / ^ 81 8(1 -y) ( (1 - z)z^ 



81 \ 1 — X z'^ \ ip' %f' 

^ f 1 _ + i ( ^(^-^) _ ^(1-^) + (-9 + %)^ + ^2 

3 V 6 V yz 

9- log(y) - 94 log(2) - 2 ( — ) log(y) log(z) 
y y^ \xyj 

^ + — ) log(x) log(y) - ( ^ + %] log(a;) log(2;) 



2 52 9(1 2(l + z)\,, ,^ , 
- + ^ - (l«g(l - ^) log(^) + Li2(x)) 

+ 2 f ^i^l (log(l - y) log(y) + U^i^y)) 



(log(l-z)log(z)+Li2(z)) , (B.IO) 



xyz 

9xz 4z 2z 

y^ xy 



Ae^{x,y,z)=-—\-—-—;^ + x --o - -3 + 



8y , 2/ 8 8y\ , ~ Sly 



81 y X x^ \ y"^ ) y (1 — x)2; 

1/ 9x2 ;^2x 82/ 8y2 _i + (8-9a;)x 2(-l + 6x)\ o 

~ « T ^ ^ 2 2 ' ^ 

b \ y X x^ z ) 

1 /46x 27x2 2(20 + 23(-2 + x)x)y\ , ,,22/ 1\ 

\Z2XZ/ \ y X J 

/2(-2 + x)x lOxy y2\ , , , 
+ ( ^ ^2 + ^ - ^2 J log(^) log(2/) 

/ 3x(-4 + xH-4z)\ 
- (1 -2 1 log(x) log(z) 

-2(^^±M!^log(y)log(.) 

-l + (8-9x)x 2-12x 3x(-4 + x + 4^) 
z"^ z y2 

X (log(l - x) log(x) + Li2(x)) 

, / 8y(l-z) , 2 + y(-10 + 9y) , 6y^^ r^^TW^^ 
+ ( ^2 ^ ^2 T j ^ ~ ^^^^-^ Li2(y)) 

/9x 3 8y 9xz 32; SyzX , , -r. / /t,-,-.n 

+ — + - + ^- — V log 1 - z) log 2; + L12 2 , B.ll 

y^ y x^ y^ y x^ / 
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2 / 16(^-y' -yz^ + {1- z)z^ + (l- z + ^'))\ 44 

81 "r 



+ i f 2 + - + ~ - -^ + ^^ + 9(l-t/)y ^ 2-6t/ 

6 \ y X z^ z 



(9-6y)z 9z2\ „ 1 / 6 27(1 

2^^ T P + o 68 - + ^ ^ 

y d \ 1 — a; yz 



/ x(9-12y) -5 + y\, 

V y V y J 

( 4(l-y)y -l + 4y\ 

- 2 (^1 -2 — + — - — ) log(y) log(2;) 

(1 - x) (-9(1 - xfx + 5(1 - x)(-l + Qx)y + 5(1 - 6x)y2) ' 



+ 2 + 



X (log(l - x) log(x) + Li2(x)) 



2-3x 9(1 -x)x 5 - 12x 8(1 - x)z Sz'^ 
X 2; x^ x^ 



X (log(l - y) log(y) + Li2(y)) 

2 - 3x 9(1 - .;•)•'' J - 12x 8(1 - x)y 8y2 
X y x^ x^ 

X (log(l-z)log(z) + Li2(z)) , 

(B.12) 

B0^ix,y,z) = 2Bs,ix,y,z) fori = 1,2,3, (B.13) 



Cei {x, y,z)=^(^l- ^ j (iTT - log(x)) , (B.14) 

Ce,{x,y,z) =^ (^1 - (z7r-log(x)) , (B.15) 

Ces (x, y, 2) = 3 - log(x)) . (B.16) 

The complete results, including the other helicity configurations, can again be found in the 
attachments of the arXiv version of this paper. 
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